Quantitative assessment of pinning forces and the superconducting gap in NbN thin films 
from complementary magnetic force microscopy and transport measurements 
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Epitaxial niobium-nitride thin films with a critical temperature of = 16K and a thickness of 100 nm 
were fabricated on MgO (100) substrates by pulsed laser deposition. Low- temperature magnetic force 
microscopy (MFM) images of the supercurrent vortices were measured after field cooHng in a magnetic field 
of 3 mT at various temperatures. Temperature dependence of the penetration depth has been evaluated 
by a two-dimensional fitting of the vortex profiles in the monopole-monopole model. Its subsequent fit 
to a single s-wave gap function results in the superconducting gap amplitude A(0) = (2.9 =b 0.4) meV = 
(2.1 lb 0.3) k^T^, in perfect agreement with previous reports. The pinning force has been independently 
estimated from local depinning of individual vortices by lateral forces exerted by the MFM tip and from 
transport measurements. A good quantitative agreement between the two techniques shows that for low 
fields, B < jUoHc2, MFM is a powerful and reliable technique to probe the local variations of the pinning 
landscape. We also demonstrate that the monopole model can be successfully applied even for thin films 
with a thickness comparable to the penetration depth. 

PACS numbers: 74.70.-b 74.78.-w 74.25.Wx 68.37.Rt 



I. INTRODUCTION 

Vortex pinning is an important characteristic of type II su- 
perconductors that allov\^s tuning its properties on purpose 
w^ithout changing the chemical composition. Hence, the in- 
terpretation of pinning mechanisms/ the search for artificial 
defects v\^ith high pinning potentials,^ and commensurable 
pinning effects by ordered arrays of defects^"^ remain in the 
focus of basic research and application-based engineering. 
On the one hand, many high-pow^er applications require 
materials with high pinning,^ i.e. high critical current den- 
sity. On the other hand, logical applications (i.e. fluxonics 
devices)^ benefit from low-pinning materials, which can be 
locally modified by introducing strong pinning sites in order 
to tune the dynamics of vortices and rectify their motion. 

A common way to investigate the pinning strength is to 
measure the sample's response to an applied magnetic field 
or current, using magnetometry or transport measurements, 
respectively. These global methods probe the average value 
of the pinning force in a material including the collective 
dynamics of the elastic vortex lattice. However, the local 
modulations of the pinning landscape originating from dif- 
ferent natural or artificial defects remain inaccessible to 
these techniques. This challenge can be addressed by lo- 
cal imaging methods, such as low-temperature magnetic 
force microscopy (LT-MFM), which is capable of correlating 
the superconducting (SC) vortex positions with the distri- 
bution of micro- or nanostructural defects. It effectively 
combines the non-invasive imaging of flux lines with the 
ability to manipulate individual vortices by the stray field of 
the magnetic tip, offering a direct access to the local pinning 
force. Nevertheless, reconciling the results of local and 
global measurements often represents a challenge. 

In this paper, we will estimate the pinning forces in nio- 
bium nitride (NbN) thin films using two complementary 
methods. Local measurements by LT-MFM will be directly 
compared to transport measurements. Thin films of NbN 
have been chosen due to their high critical temperature. 



^ 16-17 K, that makes them suitable for LT-MFM studies 
in a wide temperature range. Moreover, this conventional 
superconductor has attracted attention due to its recent 
applications in sensitive SC bolometers. 

We will also evaluate the temperature dependence of the 
magnetic penetration depth, A(r), by performing a two- 
dimensional fit of the vortex profiles within the monopole- 
monopole model. Subsequently, we will use these values to 
estimate the superconducting energy gap of the NbN film, 
which we will compare to the direct tunneling measure- 
ments from the published literature. 



II. EXPERIMENTAL DETAILS 

Epitaxial NbN thin films with a thickness d = 100 nm 
were fabricated on single-crystalline MgO (100) substrates 
by pulsed laser deposition (PLD) using a Nb (99.95 %) tar- 
get in N2 atmosphere with a pressure of 5 • 10~^ mbar. The 
base pressure in the chamber was 10~^ mbar. Prior to the 
film preparation, the Nb deposition rate was measured with 
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Fig. 1 . (a) Zero-field resistive SC transition with a full width of 
^ 0.25 K. (b) Second critical field, iXqH^2^ as a function of the 
reduced temperature, t = T/T^. Points are the experimental values 
determined from transport measurements, whereas the solid lines 
are empirical fits.^^ 



an Inficon XTM/2 rate monitor. We used a KrF excimer laser 
{Lambda-Physik) with a wavelength of 248 nm and a pulse 
duration of 25 ns. The substrates were heated up to 500°C 
during deposition. As the on-axis PLD process leads to the 
formation of droplets on the surface, which pose a severe 
problem to the MFM scanning tip, a polishing technique 
was applied to remove these obstacles, providing a peak- 
to-valley roughness below 5 nm. X-ray diffraction patterns 
showed (000 peaks (simple cubic structure) similar to previ- 
ous reports. The best samples exhibit a (offset) at 16 K 
and a sharp resistive SC transition with a width aT^ ^ 0.25 K 
[Fig. 1 (a)]. The temperature dependence of the second crit- 
ical field, /ioHc2, determined from transport measurements, 
is shown in panel (b) of the same figure. 

The LT-MFM measurements have been performed using 
a commercial scanning-probe microscope {Omicron Cryo- 
genic SFM)}^ We have used an MFM cantilever (Nanoworld 
MFMR) that possesses a force constant A: 2.8N/m and a 
resonance frequency /q ^ 80 kHz. For the transport mea- 
surements, a lOOjUm-wide bridge was structured by optical 
lithography and ion-beam etching. Transport measurements 
were performed in the standard four-point configuration 
using a 9 T physical property measurement system (PPMS) 
by Quantum Design. 

III. MONOPOLE-MONOPOLE MODEL 

The magnetic moment of the MFM tip (see sketch in 
Fig. 2) can be well approximated to the first order by a 
magnetic monopole characterized by a "magnetic charge" 
m located at a distance 5 from the sharp end of the tip 
pyramid. The field distribution from a single flux quantum, 
00 = 2.07 X 10~^^Tm^, measured just above the surface 
of the superconductor, is also similar to the magnetic field 
emanated by a magnetic monopole of 2(^o^ located at the 
depth Agff = 1.27A below the surface, with A being the 
magnetic penetration depth. Hence, the magnetic induction 
of the vortex, B(r,z), taken at a distance z above the surface, 
can be approximated by 



B(r,^): 



00 (r-ro) + (2 + 1.27A)e, 
2^ ((r-ro)2 + (2 + 1.27A)2) 



3/2' 



(1) 



where is the unit vector orthogonal to the film, Yq is the 
position of the vortex core and r is the radial distance from 
its center. This leads us to the tip-vortex interaction force 
F(r,^) in the monopole-monopole model. 



F(r,2) = mB(r,2). 



(2) 



Taking into account that the shift of the resonance frequency 
of the cantilever measured by MFM, a/, is proportional to 
the z -derivative of the normal component of the force that 
acts between the tip and the sample,^"^'^^ d^F^{Y,z), one 
obtains the following expression for the measured signal: 

fo Mo (r-ro)^-2(^ + 1.27A + 5)^ 
^•^ 2k 2n ((r -ro)2 + (2 + 1.27A + 5)2)5/2- 

As the magnetic induction of the vortex is maximal at the 
center and decays rapidly with r, the strongest interaction 
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Fig. 2 (color online). A schematic of the MFM imaging procedure 
and an illustration of the monopole-monopole model. The magne- 
tized MFM tip, driven by a piezo element, scans above the surface 
of the sample at a given distance z. During measurement, the 
feedback loop, which is typically used to stabilize the resonance 
frequency of the tip during topographic imaging, is deactivated. 
One measures a shift of the resonance frequency, a/, induced 
by the magnetic field of the vortices, B(r,2). In the monopole- 
monopole model, described in the text, both the tip and the vortex 
are approximated by magnetic monopoles at distances z + 5 and 
Agff from the surface of the sample, respectively. 



in z direction between the tip and the vortex occurs when 
the tip passes the center of the vortex. Thus, the maximal 
z-component of this force is reached at r = ToI 



max(Fj,^ 



m0o 



271 (2 + 1.27A + 5)2 



(4) 



The benefit of the monopole-monopole model is that all 
spatial parameters of the problem {z, Agff and 5) enter Eqs. 3 
and 4 additively hence the sum in the denominator can be 
redefined as an effective tip-sample distance w = z-\- A^ff + 
5. It represents the distance between imaginary magnetic 
monopoles within the tip and the vortex, as illustrated by 
white dots in Fig. 2. 

The ratio between the lateral and vertical forces that act 
during scanning typically varies from 0.3 for a tip with a 
less sharp pyramid^^ to 2/3 a/3 ^ 0.38 for the monopole- 
monopole model. ^^'^^ This results in the following maximum 
lateral component of the tip-vortex interaction force: 



max(Fiat) ^ 0.38 • max(Fj. 



(5) 
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(a)T=8.4K = 50%T, 




(b)T=10.0K = 62%T, 



(c) 7= 11.8 K = 75% T, (d) 7= 13.7 K = 87% 7, 




Fig. 3 (color online) . Vortex images measured at a distance z = 
30 nm from the surface of a NbN film after field-cooling in —3 mT at 
various temperatures: (a) 50% T,, (b) 62% T,, (c) 75% T,, and (d) 
87% T^. The slow scanning direction is top to bottom. Depinning 
of vortices by the magnetic tip can be seen in panels (c) and (d) . 



Obviously, non-invasive imaging of vortices by MFM is pos- 
sible only as long as the vortices are pinned. The tip-vortex 
interaction force can be accurately tuned during scanning by 
varying the tip-sample separation z}^ If this force exceeds 
the pinning force of an individual vortex at a natural or 
artificial defect, the vortex can be dragged aw^ay from its 
initial position. Likewise, if the tip-sample distance is kept 
constant, an increase in temperature can lead to the local 
depinning of individual flux lines during scanning due to 
the temperature dependence of the pinning force. 

IV LOCAL DEPINNING OF INDIVIDUAL FLUX LINES 

We have imaged the temperature dependence of the vor- 
tex distribution after field-cooling the sample to ~ 8 K in a 
vertical applied magnetic field of /io^z = — 3mT (Fig. 3), 
and subsequently increasing the temperature in steps of 2 K. 
The tip-sample distance w^as kept constant at z = 30 nm, 
so that Eq. (4) can be applied. The negative direction of 
the field corresponds to the magnetic repulsion between the 
MFM tip and the vortex, and therefore the vortices appear as 
red (dark-grey) objects in the false-color images presented 
in Fig. 3. In all panels, the white circles depict the positions 
of the vortex cores at base temperature, to emphasize the 
tip-induced changes in their position as the temperature 
is increased. These positions were determined by a two- 
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Fig. 4 (color online) . Temperature dependence of the averaged 
vortex profile. Blue points are the measured signal, solid lines are 
least-squares fits to the monopole-monopole model. 



dimensional (2D) fitting procedure^^ applied to the regions 
indicated by rectangles. The benefit of this method is that, 
in principle, it allows for a subpixel resolution of the fitting, 
given that the noise in the measured data is sufficiently low. 

At temperatures not exceeding ~ 50% [Fig. 3 (a)], non- 
invasive imaging of vortices takes place, indicating that the 
pinning exceeds the lateral thrust of the MFM tip. Statistical 
image analysis, such as described by Inosov et al,^^ reveals 
that the vortices form a highly disordered hexagonal lattice 
due to the pinning by natural defects. 

At higher temperatures, the decreasing contrast of the 
vortex profile signifies a natural increase in the penetration 
depth. A, that characterizes the decay of the magnetic field 
outside of the vortex core. At 62% [Fig. 3(b)], most 
vortices are still in their original positions, implying that 
the tip -vortex interaction force is still lower than the typical 
pinning force of a single vortex. Only 2 out of 14 vortices, 
visualized in the figure, have been irreversibly dragged away 
from their initial positions (white circles) to the nearest 
pinning sites with higher pinning potentials. This indicates 
the existence of a slightly modulated pinning landscape in 
the NbN film and, hence, a spatial variation of the pinning 
force, as expected for natural defects. 

At 75% [Fig. 3 (c)], the movement of nearly every vor- 
tex by the MFM tip can be seen. Indeed, because most 
vortices are irreversibly dragged away by the tip as it passes 
close to the core, such vortices appear half-cut in the im- 
age. Consequently, at this temperature the pinning force for 
the majority of the vortices is equal to the maximal lateral 
force exerted by the MFM tip onto the vortex, Eq. (5) . Only 
one vortex at the bottom-right part of the image remains 



Temperature (K) : 


8.4 


9.1 10.0 


11.8 


w (nm) 


305(11) 


308(10) 316(10) 


383(23) 


A/(r = ro) (Hz) 


0.182(4) 


0.134(3) 


0.138(6) 


max(Fi,J (pN) 


0.74(3) 


0.56(2) 


0.70(5) 



Table I. Temperature dependence of the effective tip-sample dis- 
tance in the monopole-monopole model, w = z -\- 1.27X-\- 5; the 
fitted peak amplitude from Fig. 4, A/(r = Tq); and the correspond- 
ing lateral tip-vortex interaction, max(FiaJ. 
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Fig. 5 (color online). Temperature dependence of A(r) + 5/1.27 
values (plotted on the right scale) in comparison to A(0) (left scale). 
The solid line is a fit to an empirical model with a single s-wave 
gap.^^ The value of A(0) for our film thickness was obtained by an 
interpolation of the d -dependent literature data,^^'^"^ as shown in 
the inset. The sohd grey Hne is an empirical fit. A good agreement 
of our results with the published value indicates that the 5/1.27 
correction is negligible in our case. 



stable, evidencing the locally enhanced pinning force at 
this position. On the other hand, three other vortices are 
fully dragged aw^ay as soon as the tip starts crossing their 
field lines. At even higher temperatures [Fig. 3 (d)], vor- 
tices can no longer be imaged. Here the vortices are being 
continuously dragged by the tip during scanning. 

For a quantitative analysis of the vortex profiles, their core 
positions Tq were first determined using 2D fitting. To gather 
sufficient statistics for the application of the monopole- 
monopole model, the signal a/ from w^ithin the neighbor- 
hood of every vortex (w^hite rectangles in Fig. 3) has been 
plotted vs. |r — TqI, as show^n in Fig. 4. In this figure, every 
data point corresponds to a pixel in the original MFM im- 
age, whereas data points originating from different vortices 
are combined in one plot. The resulting clusters of points 
can be fitted to Eq. (3) (solid lines) to obtain the average 
value of w = z + 1.27A + 5 for every temperature with a 
sufficiently small statistical error. The results of these fits 
are summarized in Table I. 



V LOCAL PINNING FORCES AND THE PENETRATION DEPTH 

Now we can proceed to the quantitative estimation of the 
pinning force. Combining Eq. (3) through (5), we obtain 



kw 

max(Fi,J ^ 0.38 — • A/(r = Tq). 

Jo 



(6) 



Substituting the fitting results from Fig. 4 and the known 
parameters of the tip into this expression, one can calculate 
the temperature-dependent lateral tip-vortex forces that are 
listed in Table I, with an average value of 0.67 =b 0.09 pN.^^ 



As we already know that the local pinning force decreases 
sufficiently to allow depinning of most vortices at T ^ 
12 K = 75% Tc, we can use the calculated value as an esti- 
mate of the mean local pinning force at this temperature. 



Fp (r = 12 K) (0.67 ± 0.09) pN. 



(7) 



The results of the same fit also provide a local probe 
for the temperature-dependent penetration depth,^"^'^'^'^^'^'^ 
in our case given by A(r) ^ [w(r) — z — 5]/ 1.27. The 
first temperature-independent parameter z <: X is fixed 
during measurement and known, so it can be easily sub- 
tracted (Fig. 5). The second parameter, 5, being a prop- 
erty of the tip and dependent both on temperature and the 
stray field of the sample, usually has a large uncertainty 
and requires a special calibration of the tip in order to be 
determined.^^ It may lead to a non-negligible constant offset 
of the measured penetration depth from its true value. To 
quantify the 5/1.27 correction in the present work, we re- 
sorted to a comparison of our temperature-dependent data 
(shown by spheres in Fig. 5) to the low-temperature value of 
A(0) ^ 205 nm (dark-gray diamond in Fig. 5). The latter has 
been obtained by interpolating the directly measured A(0) 
values from the literature for different film thicknesses,^^'^"^ 
as shown in the inset of Fig. 5, and taking the intermediate 
value at d = 100 nm. Good agreement between this refer- 
ence value of A(0) and our A(r) + 5/1.27 points indicates 
that the 5/1.27 correction is negligible in our case. 

The solid line in Fig. 5 is an empirical fit to the formula of 
Evtushinsky et al. which gives an analytical relationship 
between A(r) and the SC gap A(0). For a conventional su- 
perconductor with a single isotropic s-wave gap, it becomes 
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Fig. 6 (color online). Dependence of the critical current density, 
Jc on applied magnetic field for different temperatures: 8K, lOK, 
12 K and 14 K. 
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FpViV (pN) 


2.07 1.04 


0.62 


0.42 



Table II. Critical current density and average pinning force per 
vortex at different temperatures evaluated from the transport data 
in Fig. 6 for B = 0. 



A(r) = A(0) 



1 -M 



A(r) 



-1/2 



(8) 



w^here A(0) depends only on the band structure, w^hereas 
all the temperature-dependent quasiparticle effects are in- 
cluded in the approximant function^^ 

M(t) = 4(e^/2 + e-^/2)-2 ^ ut/8 + 1/(1 + ut/8). (9) 

The temperature dependence of the SC gap in Eq. (8) is 
approximated by^^ 



A(r) = Ao tanhj^^/ryr^ 



(10) 



The resulting fit yields a value of A(0) = (2.9 =1= 0.4) meV 
= (2.1 lb 0.3) k^T^, w^hich perfectly agrees w^ith direct tun- 
neling measurements^^'^"^ and is slightly above the w^eak- 
coupling limit of lJ6k^T^ predicted by the Bardeen-Cooper- 
Schrieffer (BCS) theory.^^ 

VI. GLOBAL ESTIMATE OF THE PINNING FORCES 

While MFM provides access to the pinning force of indi- 
vidual vortices, the global characterization methods, such 
as transport or magnetization measurements, explore the 
collective behavior of the flux-line lattice. They evaluate the 
mean pinning force within the v\^hole sample volume, con- 
sidering also the elastic interaction between individual flux 
lines as well as the collective pinning. For small magnetic 
fields, B <C Mo^c2^ the distance between vortices is larger 
than A. Such vortices can be treated as independent non- 
interacting objects. In this limit, collective effects can be 
ignored and the pinning force per vortex can be calculated 
as F^/iV, where is the average pinning force and N is the 
number of vortices within the sample surface. 

The pinning force is equal to the maximal sustainable 
Lorentz force that does not move vortices while the current 
flows : 

Fg(B) = VJ^B = SdJ^B, (11) 

where is the critical current density, S is the surface area 
of the sample and d is the sample thickness. The number of 
vortices is N = BS /(j)Q, hence the pinning force per vortex is 



F^/N=J,d^O' 



(12) 



Dependence of on the applied magnetic field for tem- 
peratures between 8 and 14K is presented in Fig. 6. The 
resulting temperature dependence of the pinning force per 
vortex F^/N calculated from these J^H) curves is given 
in Table II. One can immediately appreciate the agree- 
ment between the value of F^/N = 0.62 pN that resulted 
from the transport measurements at T = 12 K with that of 
(0.67 =b 0.09) pN that we extracted earlier from the MFM 
data at a similar temperature. Taking into account that the 
pinning force varies by nearly a factor of 5 in the studied 
temperature range, such an agreement within 8% between 
local and global measurements is indeed remarkable. 



VII. SUMMARY 

To conclude, we found perfect agreement between the 
values of the pinning force per vortex, estimated from local 
depinning of individual vortices by the MFM tip and globally 
from the critical current measurements. We demonstrated 
that for low fields, B <C Mo^c2^ MFM is a powerful and 
reliable method to probe the local space variation of the 
pinning landscape. The monopole-monopole model, orig- 
inally derived for d > 4A,^^ proved to be successful even 
for thin films with a thickness comparable to the penetra- 
tion depth. With this knowledge, the quality of such very 
thin films, that are actually employed for the application in 
SC bolometers, can be perfectly analyzed using magnetic 
force microscopy. 

Finally, we used accurate 2D fitting of the vortex pro- 
files to extract the London penetration depth of the NbN 
film and the SC energy gap. The statistical errors of this 
method are small enough to ensure that the extracted gap 
A(0) = (2.9 ± 0.4) meV agrees with the directly measured 
values. ^"^'^^ Although similar methods of extracting the gap 
amplitude from the muon-spin rotation (jUSR), small-angle 
neutron scattering (SANS), microwave surface-impedance 
(MSI), or magnetization measurements of the temperature- 
dependent penetration depth are well developed,^^ their 
application to the analysis of temperature-dependent LT- 
MFM images is only becoming a standard practice. 
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